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Abstract. Recent progresses in statistical mechanics indicate the Tsallis nonextensive 
thermostatistics as the natural generalization of the standard classical and quantum 
statistics, when memory effects and long range forces are not negligible. In this frame- 
work, weakly nonextensive statistical deviations can strongly reduce the puzzling dis- 
crepancies between experimental data and theoretical previsions for solar neutrinos and 
for pion transverse-momentum correlations in Pb-Pb high-energy nuclear collisions. 



1 Introduction 

Nuclear and, in general, many-particle systems are often in regimes where con- 
cepts of statistical physics come into play. Long-range forces, memory effects, 
sizeable correlations and fluctuations of the observables provide demanding test- 
ing grounds of the equilibrium and non-equilibrium statistical mechanics theory. 

In many physical applications, the kinetic description of the system is devel- 
oped in the weak-coupling limit defined by the condition Tb <C (t\, t u ), where 
Tb, t\ and t u are the duration of a binary collision, the mean free time between 
subsequent collisions and the characteristic time of the mean-field potential, 
respectively [1]. In this case, the many-body collision process is well described 
independent binary collisions in the framework of the standard Boltzmann equa- 
tion and in terms of the 8(t — t) function for the fluctuating stochastic Langevin 
force (Markovian approximation). 

If the mean-field characteristic time becomes comparable to the duration of 
a binary collision (jb < r u ), the independent binary collision approximation is 
no longer correct and the Markovian description breaks down: memory effects 
become important. 

Tsallis [2] has recently advanced a generalization of the conventional Boltzmann- 
Gibbs thermostatistics that overcomes the inability of the conventional statistical 
mechanics to tackle those many physical problems with long-range interactions, 
long-range microscopic memory, or fractal space-time constraints. There exist 
already many applications: astrophysical self-gravitating systems [3], the solar 
neutrino problem [4] , distribution of peculiar velocities of galaxy clusters [5] , cos- 
mology [6] , many-body theory, dynamical linear response theory and variational 
methods [7]. 
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In this paper, after a brief review of the Tsallis thermostatistics in Sec. II, 
we discuss (Sec. Ill) how a weakly non-ideal stellar plasma, such as the solar 
interior, can produce an equilibrium velocity distribution that deviates from the 
standard one and how this effect can be relevant to the solar neutrino problem. 
In Sec. IV, we consider the interpretation of the pion transverse-momentum 
NA49 experimental data in high-energy Pb-Pb collisions in the framework of 
the nonextensive thermostatistics. Although the physics of the two problems is 
very different, they are both characterized by memory effects and long-range 
forces and, consequently, show nonextensive statistical behaviors. 



2 Nonextensive statistical mechanics 

The principal features of the Tsallis generalized thermostatistics is based upon 
the following two postulates [2] . 

— Given pk the probability in any W different microstates k, the entropy of a 
system is defined as 

w 

1 k=i 

where q is a fixed real parameter. The generalized entropy has the usual 
properties of positivity, equiprobability, concavity, irreversibility and, in the 
limit q — > 1, is equal to the conventional Boltzmann-Gibbs entropy S = 

— The mean value of an observable O, whose value in the microstate k is Ok, 
is defined as 

w 



fc=i 



The S q entropy is nonextensive. If A e B are two independent systems Ae B, 
i.e., the probability of composite system A + B factorizes into pa+b{ u a, u b) = 
Pa(ua) Pb(ub)), the global entropy is not the sum of the entropies of the sub- 
systems, but 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) . (3) 

The single particle distribution function is obtained by the usual procedure 
of maximizing the Tsallis entropy with the constraints that the average internal 
energy and the average number of particles remain constant: 



!-(!-?) 



mv 2 
2kf 



1/(1-9) 



(4) 



When the entropic parameter q is smaller than 1 the distribution has an upper 
cut-off: mv 2 /2 < kT/(l — q) (the tail is depleted). The distribution correctly 
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reduces to the exponential Maxwell-Boltzmann distribution in the limit q — ► 1. 
When the parameter q is greater than 1, there is no cut-off and the (power-law) 
decay is slower than exponential (the tail is enhanced). 

These classical distribution can be generalized to the quantum case under 
some not too restrictive approximation (see Ref.[8] for details) resulting in the 
following analytical expression for the mean occupational number: 

{n)q = [f + (< 7 -l) / 3(£-/i)] 1 /( 9 -i)±i ' (5) 

where /3 — 1/kT, and the + (— ) sign applies to fermions (bosons), fn the limit 
q — ► 1 (extensive statistics), we recover the conventional Fermi-Dirac and Bose- 
Einstcin distribution. 



3 Non-ideal stellar plasma and solar neutrino fluxes 
3.1 Plasma parameter and solar interior 

The stellar core is usually described as an ideal plasma in the Debye-Hiickel 
mean-field approximation. However, density and temperature conditions in stel- 
lar plasmas, such as the core of the Sun, of brown dwarfs or of Jupiter, suggest 
that the microscopic diffusion of electrons and ions could be nonstandard and 
that the Debye-Hiickel approximation is not sufficiently accurate. 

The appropriate theoretical approach and the effective interactions used to 
describe a plasma can be deduced from its plasma parameter r 

r {Ze)2 ™ 

where a — n -1 / 3 is of the order of the interparticlc average distance (n is the 
average density) . The plasma parameter is a measure of the ratio of the mean 
(Coulomb) potential energy and the mean kinetic (thermal) energy. On the basis 
of its value, we can distinguish three different regimes. 

— r <C 1. The plasma is described by the Debye-Hiickel mean-field theory as 
a dilute weakly- interacting gas. The screening Debye length 



/ kT 

is much greater than the average interparticle distance a, hence there is a 
large number of particles in the Debye sphere (Nd = (47r /i)R? D ). Collective 
degrees of freedom are present (plasma waves) , but they are weakly coupled 
to the individual degrees of freedom (ions and electrons) and, therefore, 
do not affect their distribution. Binary collisions through screened forces 
produce the standard velocity distribution. 
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— r w 0.1. The mean Coulomb energy potential is not much smaller than the 
thermal kinetic energy and the screening length Rr> « a. It is not possible to 
clearly separate individual and collective degrees of freedom. The presence 
of at least two different scales of energies of the same rough size produces 
deviations from the standard statistics which describe the system in terms 
of a single scale, kT. 

— r > 1. This is a high-density/low-temperature plasma; the Coulomb inter- 
action and quantum effects start to dominate and determine the structure 
of the system. 

In the solar interior the plasma parameter is r & « 0.1; therefore, the solar 
core is a weakly non-ideal plasma where the Debye-Hiickel conditions are only 
approximately verified. 

The reaction time necessary to build up screening after a hard collision can be 
estimated from the inverse solar plasma frequency t p i — lu^ 1 = y/m/ (Anne 2 ) w 
10 -17 sec, and it is comparable to the collision time t co u — (crvn)~ l « 10~ 17 
sec. Therefore, several collisions are likely necessary before the particle looses 
memory of the initial state and the scattering process can not be considered 
Markovian. In addition to many-body collisional effects, electric microficlds are 
present and the distribution and the fluctuations of these microfields must be 
carefully considered, since they modify the usual Boltzmann kinetics. 

At the light of the above considerations, we conclude that density and temper- 
ature conditions of the solar core suggest that non-Markovian memory effects and 
long-range forces are present. Because the solar core lies in a intermediate region 
of the plasma parameter, the corrections to the standard Maxwell-Boltzmann 
distribution should be small and affect mainly the more energetic particles. How- 
ever we will see in next subsection that this small deviation strongly affects the 
nuclear reaction rates and, consequently, the determination of the solar neutrino 
fluxes. 

3.2 Thermonuclear reaction rate and solar neutrino fluxes 

The nuclear reaction rates in the stellar interiors play a crucial role in the un- 
derstanding of the structure and evolution of stars. 
The reaction rate per particle pair is defined as 



where the particles distribution function f(v) is a local function of the temper- 
ature. 

The kinetic energy of the ions in the solar core is essentially the thermal 
energy /cT Q ss 1.36 keV which is far below the Coulomb barrier, therefore only a 
small number of particles in the high-energy tail of the distribution has a chance 
of reacting: this high-energy tail plays a crucial role for the reaction rates, which 
becomes very sensitive to small changes of the distribution. 




(8) 
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For small deviations from the standard extensive equilibrium statistics, we 
can write the generalized Tsallis distribution (4) as follows [4] 

f(E) ~ (fcT)- 3 / 2 e -E/kT-5(E/kT)i ) (g) 

where S = (1 — q)/2. 

The presence of even a tiny deviation from MB in the solar core produces 
large changes of the subbarrier nuclear reaction rates and, consequently, of the 
predicted neutrino fluxes. Following the general homology relationships for the 
variations of physical inputs, see for instance Ref. [9], we can estimate the effect 
of the non-Maxwellian distribution on the fluxes [4] : 

R j = ^=e-W, (10) 

3 

for the fluxes j = 7 Bc, 8 B, 13 N and 15 0, while we use the solar luminosity 
constraint to determine the pp flux, R pp = 1 + 0.087 x (1 - R Be ) + 0.010 x (1 - 
R N ) + 0.009 x (1 - R ), and keep fixed the ratio £ = <P pep /<P pp = 2.36 x 10~ 3 . 
The power indices (3j, that appear in Eq. (10), depend on the nuclear reaction 
considered and their values have been taken from Ref. [9]. In principle, there 
could be a different parameter 5j for each reaction and it should be possible to 
calculated them from the specific interactions in the solar plasma core. For the 
purpose of estimating the possible effects of this mechanism on the solar neutrino 
fluxes, we used two simple models where 5j were used as free parameters. 

The first model uses the same 8 for all the reactions; the best fit to the ex- 
perimental data give 5 = 0.005 with a corresponding \ 2 — 35. The second model 
fits two different 5's, one (5(i7)) for the p + 7 Be reaction and the other (5(34)) 
for the 3 He + 4 He reaction. The best result gives x 2 = 20 with 5( 17 ) = —0.018 
and 5(34) = 0.030. The negative value of 5n 7 ) means that the corresponding 
distribution has an enhanced tail and that the p + 7 Be reaction rate increases. 

In spite of the fact that the values of x 2 are much smaller than the ones in 
the SSM (xssm > 7 4)j they are still large: this mechanism cannot solve the solar 
neutrino problem. However, even if S is small, it has non-trivial consequences on 
the neutrino fluxes: the boron (beryllium) flux can change of as much as 50% 
(30%). 

4 Fluctuations and correlations in high-energy nuclear 
collisions 

The recent event- by- event analysis of central Pb+Pb collisions at 158 A GcV 
(NA49 coll.) has spurred great interest because of the strong suppression of the 
pion transverse- momentum fluctuations in Pb+Pb respect to p+p collisions. A 
reduction of the fluctuations appears reasonable because correlations should be 
washed out by meson-baryon and meson- meson rescatterings; however, the size 
of this suppression is not understood. 
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Gazdzicki and Mrowczyriski [10] have introduced a definition of the event- 
by-event transverse-momentum fluctuations that is independent of the particle 
multiplicity: 

where z p± = p± — p± is a single-particle measure of the deviation of p± from 
its average value, and where Z p± = X^iCP-U ~ PT) i s the corresponding total 
contribution (N is the number of the particles in the event). 

In a very recent paper, Mrowczynski [11] has calculated this correlation mea- 
sure <P P± for a pion gas in global equilibrium within the standard extensive ther- 
modynamics. His result (<P P± = 6.5 MeV at T = 120 McV) is sensibly greater 
than the experimental value <P p± = 0.7 ± 0.5 MeV, which has been measured in 
the central Pb+Pb collision by the NA49 collaboration [12]). 

What is the origin of such a large discrepancy between theoretical results 
and experimental data? Can the conditions of density and temperature at the 
early stage collisions modify the correlations between the produced particles? 

The quoted theoretical result [11] has been obtained within the framework 
of the conventional quantum statistical mechanics. Therefore, pions, which are 
sufficiently light to show quantum degeneracy, follow the standard Bose- Einstein 
distribution. However, it is common opinion that, because of the extreme condi- 
tions of density and temperature in ultrarclativistic heavy ion collisions, mem- 
ory effects and long-range color interactions affect the thermalization process 
[13, 14, 15]. In fact, if the quark-gluon plasma forms at the early stage of the 
collision, strong chromoelectric fields appear within the parton gas. Recent in- 
vestigations, both for for QED and QCD, have shown that these strong fields 
are connected to the presence of non-Markovian process in the kinetic equation 
[16]. 

The foregoing considerations suggest to calculate the equilibrium correlation 
measure in the framework of the nonextesive statistics. 

On the basis of the generalized thermodynamics relations [2] , it can be shown 
that 



z 2 

^ p 



where (n) q is the mean occupation number of Eq.(5) and 
(Z 2 P1 ) 1 f d* P 



(N) 



4/(^(^) 2( ^' (13) 



where 



and 



1 f d 3 p f d 3 p 

^ = - p JWf> pl <">' • " = J (2iF < " ) ' ( ' 
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are the particle fluctuations in nonextensive (q ^ 1) statistics. For q = 1, we 
recover the well known fluctuations expression for fermions (— ) and bosons (+). 
From Eq. (15), we see that the expression for the generalized fluctuations has the 
same structure of the standard one modified by the factor 1/[1 + (1 — q)f3(E — /i)]. 

For the boson case a value q > 1 implies a very strong suppression of <1> P± . 
If we consider a pion gas and fix the freeze-out temperature to T = 120 MeV 
(obtained from the analysis of single particle spectra), we reproduce the experi- 
mental value {$ p± = 0.7 ± 0.5 MeV) using q = 1.015. Hence, a small deviation 
from the standard statistics (q — 1 = 0.015) is sufficient for eliminating the 
puzzling discrepancy between theoretical calculations and experimental data. 

Finally, we observe that fluctuations are the more strongly modified the larger 
the mass of the detected particle for a given value of the deformation parameter 
q. Hence, this model predicts stronger nonextensive effects for correlations of 
heavier mesons and baryons: future measurements of <P p± for these particles 
should find an even larger reduction when plasma is formed. New data and 
investigations are necessary to gain a deeper understanding of the high-energy 
heavy-ion observables. 

5 Conclusion 

We have considered the relevance of the generalized nonextensive statistics to 
the the solar neutrino problem and to the interpretation of the pion transverse- 
momentum correlations in Pb-Pb high energy collisions experimental data. Al- 
though the two problems involve rather different physics, both appear to show 
memory effects and long range forces. Discrepancies between experimental data 
and theoretical models can be strongly reduced when weakly-nonextensive (\q — 
1| < 0.02) statistical effects are considered. 

In the solar neutrino problem a small deviation from the Maxwcll-Boltzmann 
distribution produces strong modifications of the thermonuclear reaction rates 
and, consequently, modifies neutrino fluxes by amounts comparable to those that 
constitute the solar neutrino problem without affecting bulk properties such as 
the sound speed or hydrostatic equilibrium. 

In high-energy ion collisions, nonextensive statistics predicts a reduction of 
the correlation measure <P p± in Pb-Pb collisions ($ p± = 0.7 ± 0.5 MeV) respect 
to p-p collisions (<P p± = 4.2 ± 0.5 MeV) due to plasma effects. This reduction 
has been seen in the experiments and could be interpreted as a signature of the 
transition to a quark-gluon-plasma phase driven by the extreme conditions of 
density and temperature in the early stage of the collision. 
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